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We study a model based on the dihedral group D4 in which the dark matter is stabilized by the 
interplay between a remnant Z2 symmetry, of the same spontaneously broken non-abelian group, 
and an auxiliary Z| introduced to eliminate unwanted couplings in the scalar potential. In the 
■») ■ lepton sector the model is compatible with normal hierarchy only and predicts a vanishing reactor 

£S) \ mixing angle, 813 — 0. Since m vi = 0, we also have a simple prediction for the effective mass in 

■ terms of the solar angle: = \m U2 \ sin 2 #0 ~ 10~ 3 eV. There also exists a large portion of the 

model parameter space where the upper bounds on lepton flavor violating processes are not violated. 
We incorporate quarks in the same scheme finding that a description of the CKM mixing matrix is 
possible and that semileptonic K and D decays mediated by flavor changing neutral currents are 
under control. 
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I. INTRODUCTION 



We have strong evidence about the existence of dark matter (DM) [HQ. A good DM candidate must be neutral and 



> 

stable or with a decay length bigger than the age of the universe and give the correct relic abundance [3|. There are 

r — ■ . 

several extensions of the standard model predicting good DM candidates; however, it turns out that in many models 
the stability of the DM is obtained introducing ad- hoc assumptions, see for example the review Q|. Any of these 
models may be correct but certainly it would be desirable to provide a fundamental explanation of the origin of the 
stability. In Q it has been pointed out that the stability can be guaranteed by a residual Z2 symmetry arising from 
the spontaneous breaking of a non-abelian flavor symmetry; the same Z2 also acts in the neutrino sector and has a 
strong impact on the phenomenology of neutrino masses and mixing. In that model the flavor symmetry is the group 
of the even permutations of four objects A4 whose irreducible representations are three singlets and one triplet. To 
avoid a direct couplings to quarks and charged leptons, the DM candidate is assigned to a triplet representation, while 
leptons and quarks to singlets of A4. After electroweak symmetry breaking, A4 is broken into its subgroup Z2 under 
which two component of the triplet DM are automatically charged; eventually, this prevents dangerous couplings with 
the Higgs fields of the model. Such an idea has been then further studied and extended in refs. 0,0). 

The interplay between decaying dark matter and non-abelian discrete flavor symmetries has been considered in a 
number of subsequent papers; for instance, in [8-10| non-abelian discrete symmetries prohibit operators that may 
induce too fast dark matter decay; in ll| a non-abelian discrete symmetry (not a flavor symmetry) has been used 
to stabilize the scalar DM candidate (similar to what has been discussed in the inert scalar models [l3|) and the 
matter sector has not been considered. Therefore the models in and [ll| are substantially different to the idea 
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introduced in 




In this paper we adopt the point of view elucidated in [5-7], studying a flavor model where the stability of the 
DM is caused by the interplay between a remnant Z 2 symmetry of the D4 group and an auxiliary which allows 
to eliminate dangerous couplings in the scalar potential. Since D4 contains only singlets and doublets it is highly 
non-trivial to still be able to generate the mechanism for dark matter stabilization; in addition, the same non-abelian 
symmetry also acts on leptons and quarks, giving acceptable phenomenology in both sectors. 

The relevant differences of our model compared to [5j-[7( can be summarized as follows: 

• for the first time, we extend such a mechanism to incorporate quarks transforming under non-trivial represen- 
tation of D4; in this framework, we are able to reproduce the correct order of magnitude of the quark mixing 
angles, a quite remarkable result; 

• charged leptons are non-diagonal (with hierarchy among the eigenvalues naturally reproduced with 0(1) Yukawa 
couplings) and completely responsible for the atmospheric mixing angle in the neutrino sector, instead of being 
diagonal as in 0-0 ; 

• although the Higgs sector is extended with three more scalar doublets and one singlet, the neutrino sector 
contains only two right-handed neutrinos. The model can be considered minimal in this respect. 

The paper is organized as follows: in section |n] we present the relevant features of the model, discussing the group 
properties of D4 and the assignments of leptons and Higgs fields to the irreducible representations of the group. In 
section IIIII we discuss the scalar potential of the theory and describe in details how the DM stability arises in our 
model; sections IIVI and [V] are devoted to the neutrino phenomenology and to the estimate of some relevant lepton 
flavor violating processes, respectively. In section I VII we discuss the quark sector and give an order of magnitude 
estimate of some of the flavor changing neutral current processes; eventually, in section fVIII we draw our conclusions. 



II. MODEL 



We assign the fields of the model into irreducible representation of D4, the dihedral group of order four 13j, see 



also [14 1 . It has five irreducible representations, four singlets li,2,3,4 and one doublet 2. The generators of the group 
fulfill the relations: 

A A = B 2 = 1, ABA = B. (1) 

The one-dimensional representations are characterized by A = B = 1 for 1 4) A = 1, B = — 1 for I2, A = — 1, B = 1 
for I3 and A = —1, B = —1 for I4. The generators for the two-dimensional representations are 
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An interesting feature of D4 is that the product of two doublets contains only singlets: given (01,02) ~ 2 and 
(&i, 62) ~ 2 we have: 

ai&2 + a 2 &i ~ li, aib 2 ~ a 2 bi ~ l 2 , ^ 
ai&i + a 2 &2 ~ I3, CLih - a 2 b 2 ~ U- 

For the singlets: lj x lj = li, li x lj = lj for i = 1, • • • , 4, I2 x I3 = I4, 1 2 x 1 4 = I3 and I3 x 1 4 = 1 2 . The standard 
model Higgs doublet is taken as a singlet 1 4 ; we assume three further Higgs doublets, one of them transforming as 
a singlet I3 (H') and the other as a doublet of D4, r\ = (^1,^2) ~ 2. In order to correctly describe both lepton 
and quark sectors, we need to introduce a scalar SU(2) singlet flavon (f> in the I2 representation. Two right-handed 
neutrinos (N\, N 2 ) in the doublet representation Nn are necessary ingredients to give mass to the neutrinos via the 
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TABLE I: Assignment of the lepton and Higgs fields under SU(2), Da and Z 2 - 



type-I seesaw mechanism. The assigment to the irreducible representations of D4 as well as the charges under are 
listed in Tabffl 

The invariant D4 x Z( Lagrangian in the lepton sector is as follow: 

ylL e (N D rf) Xl + y%L^N D r,) l2 + y» 3 L T (N D v)i 3 + (4) 
+M 1 N D N D + h.c. , 

where A is a large energy scale. Here we have only considered terms up to one-flavon insertion. It turns out that 
higher powers of Higgs fields and flavon insertions can only modify (to a negligible level) the couplings y\ but are not 
able to generate new Yukawa interactions. In the neutrino sector, due to the Z^ symmetry, the previous lagrangian 
is modified by operators with at least two powers of <p, that we assume here negligible. We will investigate in details 
its phenomenological consequences in sections IIVI and |Vl 

III. SCALAR SPECTRUM AND STABILITY OF THE DARK MATTER CANDIDATE 

The invariant scalar potential is of the form: 

V = fJ? H H*H + n 2 H ,H*H' + ulirfrfi^ + /i<^ 2 + 

+ M'7 t »?)ii(»7 t *?)i 1 +*2(»? t r?)i. 2 (i7 t ?7)i 3 + A 3 (?7 t 77)13 (77+77)13 + A 4 (?;t 77 ) l4 ( ?7 t, ? ) l4 + 

+ A'lfaVHiMii +A' 2 (?7V)i 2 (w)i 2 +A' 3 (j7V)i 3 (w)i3 +A^V)i 4 (w)i4 +A 5 (r7 + r?)i 1 ff t #+ 

+ \' 5 (rfH)(Hl v ) + \ 6 [(rfrf) lt HH + h.c] + \ r (rfri) ll H , *H' + \' 7 (rfH')(H*r))+ (5) 

+ ^{{rfrf^H'H' + h.c] + \ 9 (H'^H')(WH) + X' 9 (H^H)(WH') + 

+ \ W [(H*H)(H*H) + h.c] + \ u (HlH)(H*H) + \ X i{H* H')(H* H')+ 

We assume a vev structure of the form: 

(H) = v, (H') = v', (771) = v m , (m) = v m , (</>) = V4, (6) 

where the various vevs Vi are obtained solving the coupled differential equations dV/dvi = 0. Assuming for simplicity 
real vevs, we have carefully checked that, for suitable parameter choices of the potential V, an allowed local minimum 
is: 



(7) 
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which is the crucial point to justify the stability of the DM based on symmetry arguments. After electroweak symmetry 
breaking we can write: 



h = 



v v + Vi + 




m 



H' = 



\v n + V2 + iA 2 
H'+ \ 
v' + H' + iA' j 



(8) 
(9) 



and the physical spectrum involves four neutral scalars, three pseudoscalars and three charged scalars (plus one 
Savon). To maintain the notation compact and avoid unnecessary complications, we work in the limit of decoupled <p 
(that is £1,2,3 = 0), which does not modify any of the results discussed in the paper (and, of course, the vev alignment 
(77) ~ (1, 1)). In such a limit, the mass matrices of the three sectors (S=scalar, A=pseudoscalar, ff + =charged) can 
be generically written in the following way: 



(M 



S,A,H+\2 



M^ A ' H+ 
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23 

Mli A > H+ 
Mli A - H+ 



Mli A > H+ 



(10) 



The relevant feature here is that the 2x2 sub-block corresponding to the 3 — 4 sector is symmetric and can be put 
in a block diagonal form by a maximal rotation. This corresponds to a rotation in the corresponding bidimcnsional 
subspace which defines the mass eigenstates of the subsector. After this change of basis, we are left with block-diagonal 
mass matrices made by 3 x 3 matrices (one for scalars, one for pseudoscalars and one for charged), and lxl blocks 
corresponding to the isolated DM sector: 
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(11) 



where 



f S,A,H + \2 
( m DM ) 



M S,A,H+ _ M S,AM + 



(12) 



The explicit expressions of the physical masses are complicated functions of the potential parameters and their expres- 
sions do not reveal any important features to be mentioned here beside the fact that, as expected, the pseudoscalars 
and charged mass matrices have a zero eigenvalues corresponding to the Goldstone bosons. 

We have verified that, for a suitable set of the potential parameters, m% M is the lightest mass of the neutral 
stable states and it is smaller than the corresponding charged states carrying the same Z 2 parity; in addition, all the 
masses arising from the diagonalization of the 3x3 matrices in eq. (TTT|) are positive and do not violate the electroweak 
constraints on the T, S and U oblique parameters 15]. 



We are now in the position to discuss the mechanism for the DM stability in our model. The condition in eq.© 
can be rewritten as: 



(H) 



(v) = 



(13) 



The generator B of D4 acts as the identity on the singlet representations (li for H and I3 for H 1 ). On the other 
hand, from equation ([2]) we also see that B, which is the generator of a Z 2 symmetry, leaves invariant the vector (rj). 
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In a compact 4x4 form, such a generator can be written as: 



1 
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\0 
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(14) 



where it is understood that the first two entries work on the H and H fields and the others on the components of 
the r\ field. After performing the 3 — 4 rotation that diagonalizes the corresponding entries in the mass matrices, the 
generator can be cast in a diagonal form: 

^10 \ 
10 
10 
\ -1 / 

showing the peculiar feature of the negative value in the (4, 4) entry. The previous 2x2 rotation defines the mass 
eigenstates in the D4 doublet subspace of eq.®: 



(15) 



(16) 



that, under the remnant Z 2 symmetry defined in ea.(|15|). transform as 

Vp +Vpi Vrn -Vm , 

with vevs 1 

(r/m) = 0, (rip) = V2v v . 
The same conclusions can be drawn working in a different basis where 

1 A „ / 1 



A = 



-1 



B = 



-1 



(17) 



(18) 



(19) 



In that case, r\ develops a vev along the direction (1,0) and the component that does not take vev is charged under 
the remaining Z 2 parity 2 . 

Let us now make more transparent the role of the Z%. If we introduce linear and/or trilinear terms in the scalar 
potential transforming as I2 and I4 representations (like 4> in our case), the vev of r] would not be aligned along the 
direction (1, 1) and the residual Z 2 would be broken, causing the decay of the DM. In fact, one can add contractions 
of the form (rfr])i 2 into I2 singlet of D4 and the minimizing equations for the rj components would admit solutions 
only if 7^ (772); the auxiliarly symmetry Z 2 (under which tj> — > —<fi) only allows quadratic and quartic terms in <fi 
which transform as a singlet li and avoid the dangerous Z 2 -breaking contractions. 

We see that r] m is the only Higgs field charged under the Z 2 symmetry. This prevents any coupling with other 
Higgs fields and, considering also that the original t) field does not couple to quarks and charged lepton bilinears (but 
only to heavy right-handed neutrinos) , the component of r\ m corresponding to the lightest ^2-odd neutral spin zero 
particle is the DM candidate of the model, namely the combination 77^ — r)' 2 . 

Since the DM couplings with the standard model particles are very similar to the ones in ref. Q , we expect almost 
the same DM phenomenology as described in [jj, with rriDM in the range few GeV < vtium < 100 GeV. 



1 For a similar change of basis in Qq see [lfij . 

2 We thank Luis Lavoura to point out this possibility. 
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IV. CHARGED LEPTONS AND NEUTRINOS 



In the neutrino sector, the Dirac and Majorana mass matrices derived from eq.(j4]) read: 



( v\ vl 
y v 2 -vZ 
Vj/3 vZ 



Mr = Mi 



1 

1 



(20) 



and the resulting light neutrino mass matrix, from type-I seesaw mechanism m„ = —moM R m D , can be parametrized 
as: 



2A 2 2AC 
in,, = I -2B 2 
2AC 2C 2 

where A 2 = (y^v,,) 2 /M 1: B 2 = (y^v v ) 2 /Mi and C 2 = (y%v n ) 2 /Mi. This matrix is diagonalized by: 

Uj ■ m v ■ U v = D u ; 

to find Uu, we first compute the eigenvalues \m Vi \ and eigenvectors of m v , related in the following way: 



(21) 
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(23) 



The zero eigenvalue can be assigned to m Ul or m U3 . In the latter case, for any ordering of the remaining mass 
eigenstates, the solar angle cannot be reproduced. On the other hand, for m Ul = 0, we get an appropriate U u if we 



associate |m„ 2 | = 2(\A\ 2 + \C\ 2 ) and |? 



2\B\ 2 , with the condition \B\ 2 > (\A\ 2 + \C\ 2 ) to fit the atmospheric 



mass difference. The resulting U v is then given by: 



\MC* 
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2 0\ 



A* y/\A\ 2 + \C\ 2 Cy/\A\ 2 + \C 
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(24) 



with 



tan Or. 



(25) 



Now we consider the charged sector. Since we work in the left-right basis, the mass matrices MiMj and MlMt are 
diagonalized by the unitary matrices Ug and Vi, respectively: 



rti 



Uj-M e Mj-U e = D\ 
Vl-M\Ui-Vi = D 2 , 

where Di is the diagonal charged lepton mass matrix. 

After electroweak symmetry breaking the mass matrix of the charged leptons has the form: 



(26) 
(27) 



Mi 



I y[v'e 

y l 2 vs y\v' 'e 
\ y\v' y\v 



(28) 
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where we have defined (cj)}/A = e. The charged lepton masses are given by: 



i e « y[v'e 

2 / / /? / / 

i M » e / 2 f =^P= 29 

VK) 2 + («'y 5 ) 2 



V ~ \/ ("2/ 3 ) 2 + ( u '2/ 5 ) 2 ■ 



We can easily see that the muon mass is suppressed with respect to m T by a factor of e and enhanced with respect to 
m e by roughly a factor of v' /v, which is smaller than 1 (see below). If we concentrate only on the fi — r submatrix, 
we can rewrite it as: 

m - =(:!;)• <»» 

where a, 6, c and d are products of Yukawa couplings and vevs. We can define two unitary rotations Vl and Vr 

V L —( C ° S ° e Sm9e | Vr—( cos0r sm ° R \ ( 31 ) 
I — sin 8e cos 9i J ' I— sin Or cos 8r I 

such that: 

t ( ml \ 

VjJ M^iT Vr — " 2 . (32) 
\ y 

Given the neutrino mixing matrix in eq. (f!M)) . it is evident that the atmospheric angle originates from Vl while the 
angle 9r is a free parameter. The global 3x3 charged lepton mass matrix is then diagonalized by the following 
rotations: 

Vl = ( « r"r ) . ^ = ( " " ] ■ (33) 





To get an estimate of the magnitude of the Yukawa parameters in ([28]) we proceed as follows. The vevs v and v' are 
fixed from the minimization condition of the scalar potential, then the four Yukawa couplings 2/2345 (assumed to be 
real) are determined from the atmospheric angle, the /x and r masses and the angle Or. Let us assume 0? = ir/4 (in 
good agreement with the experimental data); then, the conditions for having a fir invariant submatrix (|30p can be 
deduced using: 

MrrKr - Vr. ( f m ° ? ) Vl . (34) 

We get: 

a 2 + b 2 = c 2 + d 2 
„ „ m 2 + m 2 

c 2 + d 2 = r 2 * (35) 



db 



i 2 -ml 



In this way we can write three out of four parameters (for instance, a, b and d) in terms of the charged lepton masses, 
the atmospheric angle (supposed to be maximal here) and c. The latter is related to the Or angle by: 

2(ab + cd) 

tan 20 R = - \ 2 J . 36 
a z — b z + c z — d z 
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It is easy to show that the system of equations in (|35|) has real solutions only of the form a ~ c and b ~ ci from which 
we deduce 3 : 



y^e ps y 3 v a j/^'e. (37) 

For 1/2 ~ 2/5 the first relation implies v'/v ~ £ whereas the second one requires a moderate fine-tuning of order 
2/3/2/4 ~ e2 - 



Finally, the lepton mixing matrix is given by: 

' cq s q 

Uiep = Uj ■ U v = sm9iSQ —CqSitiOi cosdi | (38) 

\—COs6gSQ Sq COS 6g StD.Be 

where, considering the relations in cq. (|37|) . gives: 

tan0 2 3 ~4 e - (39) 
2/3 

To reproduce the correct maximal mixing in the atmospheric sector we need 2/3/2/5 ~ £ - We clearly see that our model 
predicts a vanishing 613. Since also m\ = 0, the effective mass entering the neutrinoless double beta decay assumes a 
particularly simple expression: 



\mp{i\ = \m V2 \s% = ^Am%8%, (40) 
where Anig = \m V2 1 2 — \m Vl | 2 , with numerical values in the interval: 

0.00054 eV < \m w \ < 0.0012 eV . (41) 

V. ESTIMATE OF LEPTON FLAVOR VIOLATING PROCESSES 

Since we have more than one SU (2) Higgs doublet coupled to charged leptons, our model allows for lepton flavor 
violating processes (LFV) at tree level such as r — > 3/i and r — > fiee (see [17J for other examples of renormalizable 
models based on dihedral groups). Here we are not interested in a full study of the LFV but just to show that the 
model prediction for them can be easily maintained below their upper bounds. In the interaction basis, the Yukawa 
matrices Y and Y' can be deduced from: 

( Q \ ( y{ \ 

(42) 

V y l 3 J \ y l 5 J 

where we have reabsorbed e into y[, y\ and y\. Once we rotate the fields L and l c to the mass basis, the new Yukawa 
matrices Y and Y' read: 

Y = ul-Y-Vi, Y' = uj-Y'-V e . (43) 

They are not diagonal and contain non vanishing [i~r entries, as it can easily deduced using eq. (|3"3"|) . The Higgs fields 
should also be expressed in the mass basis but we do not take this additional rotation into account since it would 
only introduce additional mixing angles as suppression factors in the branching fraction computations (we are then 
working in the case where the LFV processes are the largest allowed in our model). 
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The relative hierarchy among the two groups of parameters depend on the choice of 6r. 
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As explained above, the Yukawa couplings yi are fixed from the value of fermion masses m e , m M and m T , the vevs v 
and i>' 4 (determined by the potential parameters) and the arbitrary angle Or of the Vt unitary matrix. To get realistic 
estimates, we performed a numerical simulation with the constraints defined below eq.(|12[). It turns out that we can 
always find solutions with v > v', which implies Y' > Y, see eas. ([57|) and (|4"2"]l . In this case, the decay width for the 
r — > 3/i process 5 is approximated by: 



T(r 3/i) 



5 ( yi yi 

6 x 2 9 7r 3 m 



(44) 



Numerical examples of the branching ratio Br(r — > 3/i) are given in TabJTIl where we choose the vev v, the ratio v/v' 
and ran 1 as independent variables. We fixed the value of Or to svo.0r — 0.9277 for which the product of the Yukawas 
Y'Y' is maximal and the branching ratios are the largest possible. As we can see, there is a region of the parameter 



v (GeV) 


v/v' 


m H > (GeV) 


m H (GeV) 


rtioM (GeV) 


Br{r -> 3/i) 


224 


3.6 


140 


115 


89 


8.1 x 10 -9 


225 


3.8 


201 


98 


87 


2.5 x 10" 9 


225 


3.8 


175 


132 


60 


3.1 x 10" 9 
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1.5 x 10~ 9 
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TABLE II: Branching ratio for the process r 
3.2 • 10" s 0y. 



3/i as deduced from our model. The experimental bound is Br(r — > 3/i) < 



space where the branching ratio for the tau decay is well below the experimental upper bound Br(r — > 3/i) < 3.2 • 10 



18| . On the other hand, the first entry is very close to the upper limit, showing that a sector of the parameter space 



will be tested in the near future at the LHC. 



VI. THE QUARK SECTOR 

In this section we discuss the extension of our model to the quark sector. The assignment of the quark fields to the 
irreducible representation of D4 is listed in Tab lllll 





Qi 


Q 2 


Qs 


<A 




Q3 


SU(2) 


2 


2 


2 


1 


1 


1 




1 4 


1 2 


li 


14 


1 2 


li 


zl 


+ 




+ 


+ 




+ 



TABLE III: Quark assignments in our model. 



The Lagrangian for the down-type quarks reads as follows: 

C down = yfQiq c 1 H + y*Q 1 q c 2 H' + yiQ 2 q°H' + yiQ 2 q c 2 H + ylQ 3 q c 3 H+ 
+ i (yiQ iq c 3 H> + y$Q 2 q c 3 H + ylQ^H' + y$Q 3 q c 2 H) . 



4 In the computation of the branching ratios, we set = v' . 

5 The process r — ¥ fiee is suppressed by Y e ' e . 
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For the up- type quarks the Lagrangian has the same structure with the obvious replacements yf — > yf and (H, H 1 ) — > 
(H, H'), where H — —1T2 H\ The mass matrices are then: 



/ u,d u,d 1 

I 2/1 v y 2 ' v 

u.d 1 u.d u.a 1 ,■ 1 , ■ ■■, 

y 3 v y 4 ' v y 7 < ve | . (46) 

u.d 1 u.d 

\ y 8 v £ v£ 




With such a texture we can easily fit the quark masses and the CKM mixing angles. In particular, given that v > v' 

{ v j "u r ICTTW c u,d • 1 ,1 u.d u.d u,d u,d j u,d u.d mi 

(as discussed before eq.(|44J) ), we can fix y x 4 5 m such a way that m l ~ y Y v, m 2 ~ 2/4 n andro 3 ~ y 5 v. then 
the Cabibbo angle is given by: 

V^4 J/4/ W 

and it can be fit to its experimental value for a suitable choice of the vev ratio v' /v (with y^' d of the same order of 
J/4' ). Taking y£' 7 of the same order of y§' , we also have U„6 « ew'/w and V c b ~ £ (which fixes e ~ 0(0.04)). 
Since the Cabibbo mixing arises from both the up and down sectors, we can have s — d and c — u tree-level transitions 
mediated by Higgses. This implies that decays like K + '° —> ir + '°ll (in the down sector) or D + '° tt +: °U and 
Df — > K + ll (in the up sector) can exceed their experimental bounds. Since the coupling Y' ee is suppressed by the 
electron mass, the pairs H can be whereas the case /i ± r =F is kinematically excluded. The tree-level transitions 

b — s, b — d are suppressed by e and we only consider B + — > K + fi + fi~ and B + — > 7r + /i + /i~. In order to give an 
estimate of such processes, we work in the worst case of unity mixings and (adimensional, that is stripped of the 
meson masses) form factors. We then have: 



Y > y> I 



rp+ -> ir+ M + /0 « O IOXJ 2 . (48) 



y/ y-/ 12 



We observe that for u ~ 224 GeV (the worst point in Tab Hi)) , the couplings 2/2,3,4,5 m ea. (l4*2l should be of 0(1O -2 ) 
to reproduce the r mass and then Y w , ~ 10~ 2 . In the down sector we have ~ 2 • 10~ 2 , yf ~ 6 • 10~ 4 , t/f ~ 3 • 10~ s 
and 2/2 ~ J/4 an d therefore ~ 6 • 10 -4 . A similar reasoning in the up sector gives Y' cu ~ 7 • 10~ 3 . Using the above 
values, we computed the branching ratios for the meson decay processes in eq. (|48p . They are summarized in Tab II VI 

We clearly see that all branching ratios are below their upper limits. We have also checked that the mass difference 
in the kaon system, driven by the K° — K° oscillation, is around 10 -14 GeV, to be compared with the experimental 
value ~ 10 -12 GeV. We stress that, even if we have taken a particular point in the parameter space to make our 
estimates, the exercise can be repeated for different input values with similar conclusions. 



VII. CONCLUSIONS 

In this paper we have discussed an extension of the standard model based on the non-abelian discrete group D4. 
We introduced three more SU (2) Higgs doublets, a combination of them giving a good dark matter candidate, one 
standard model singlet and only two right-handed neutrinos, a remarkable feature if compared with the models in 
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decay 


model prediction 


experimental bounds [181 


Bi(K + - 


-> 7T + /l + /i ) 


5.5 • 10" 9 


< 8.1 ■ 1CT 8 


Bi(K° - 


» 7rVV~) 


1.2 ■ 10~ 9 


< 2.9 ■ 1CT 9 


Br(D+ - 


-> 7T + /i + /l~) 


1.3 ■ 1(T 8 


< 3.9- 1(T 6 


Br(D° - 




5.2 ■ 10~ 9 


< 1.8- 10" 4 


Bt(D+ - 


->• K+fi+fj,-) 


8.3 ■ 10~ 9 


< 3.6- 1CT 5 


Br(B+ - 




6.4 ■ 10~ 8 


< 5.2 ■ 1CT 7 


Br(B + - 




3.2 • 10" 10 


< 1.4 • 10" 6 



TABLE IV: Branching fraction estimates for some interesting processes in our model. We fixed m H i = 140 GeV, v = 224 and 



. The stability of the DM candidate is not imposed ad-hoc but directly follows from the remnant subgroup of 
the broken D^, as we explained in details. Within the same framework, we incorporated a description of the charged 
leptons and neutrinos, showing that the normal hierarchy (with m vi — 0) and a vanishing t9 13 are natural predictions 
of our model. This also allows to get a range of values for the effective mass mpp, which turns out to be in the interval 
[0.5, 1.2] • ICR 3 cV. On the other hand, the solar and atmospheric angles are free parameters that can be easily fixed 
to their corresponding experimental values. We have carefully checked that, in a large portion of the parameter space, 
the model does not conflict with the upper bounds on some lepton flavor violating processes, like r — > 3/i. Finally, 
we extended the D4 symmetry to the quark sector, showing that the correct order of magnitude for the CKM angles 
can be easily reproduced assigning the quark fields to non trivial representation of D±. Three- level flavor changing 
neutral current processes, generated by non- vanishing off diagonal Yukawa couplings, can be maintained below their 
current experimental bounds. 
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